Order parameter and vortices in the superconducting Q-phase of CeCoIn 5 
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Recently, it has been reported that the low-temperature high-magnetic field phase in CeCoIns (Q- 
phase), has spin-density wave (SDW) order that only exists within this phase. This indicates that the 
SDW order is the result of the development of pair density wave (PD W) order in the superconducting 
phase that coexists with d-wave superconductivity. Here we develop a phenomenological theory for 
these coexisting orders. This provides selection rules for the PDW order and further shows that the 
detailed structure of this order is highly constrained. We then apply our theory to the the vortex 
phase. This reveals vortex phases in which the d-wave vortex cores exhibit charge density wave 
(CDW) order and further reveals that the SDW order provides a unique probe of the vortex phase. 

PACS numbers: 



The low-temperature high-magnetic field phase in 
CeCoIns (Q-phase) has been thought to be the best ex- 
ample of a FFLO superconductor [H, 0, S I3| and has 
thus generated a tremendous interest [5j, |f|. However, 
the recent measurements of Kenzelmann et al. [?:], sug- 
gest that this point of view should be altered. This im- 
portant experimental discovery shows that the Q-phase 
reveals itself through the appearance of an incommensu- 
rate spin density wave (SDW) order. What is striking 
about this SDW order is that it vanishes when supercon- 
ductivity vanishes at high magnetic fields. This implies 
that superconducting order is the primary order param- 
eter with the SDW order induced as a secondary order 
parameter. A possibility for such superconducting order, 
as pointed out by Kenzelmann et al. is pair den- 
sity wave (PDW) superconductivity. If the SDW order 
is associated with a wavevector Q, then the PDW order 
must have the wavevector — Q to be able to induce the 
SDW order. The SDW order has Q = (q,q, 0.5), which 
is too large to be a consequence of the long-wavelength 
modulation of a FFLO phase [1 II Q. The PDW or- 
der is more akin to the 7r-triplet staggered pairing sug- 
gested by Aperis @ or to the the PDW order suggested 
in La2- x Ba x Cu04 at x = 1/8 [13] ■ The ensuing physical 
picture is then a d-wave superconductor at low fields with 
PDW order appearing through a second order phase tran- 
sition at high fields. These two types of superconducting 
order will coexist in the Q-phase. 

The observation of this PDW order raises a series of 
deep questions about the origin of this phase. To help 
address these, we have developed a phenomenological 
theory for this PDW order. Our approach complements 
that given by Kenzelmann et al. and is based on 
irreducible representations of the full space group. We 
find that this theory strongly constrains the PDW order 
and provides useful information about the vortex phase. 

PDW superconducting order parameter: Our approach 
is to classify the PDW order in terms of irreducible repre- 



FIG. 1: Directions of Q i used in the text. The field is applied 
along the direction Q 1 



sentations of the full space group [l3| • For CeCoIn5 this 
is P4/mmm. For order appearing at a wavevector Q, 
the order parameter is defined by the irreducible repre- 
sentations of Gq (set of elements conserving Q) and the 
star of the wavevector Q (set of wavevectors symmetri- 
cally equivalent to Q). For Q — (g, q, 0.5) = (q,q, —0.5) 
Gq = {E, C2 n ,cr z ,o-Q} with the 180°-rotation around 
the axis (1,1,0), o~ z and o~q the mirror operations at the 
basal plane and the plane perpendicular to (1,-1,0), 
respectively. Note (0, 0, 1) is a reciprocal lattice vec- 
tor. In Table I, we give the irreducible representations 
of Gq together with representative basis functions for 
spin-singlet pairing (scalar functions ip(k) [ill]), spin- 
triplet pairing (vector functions d(k) [U|), and spin 
density order (Si). To define the additional order pa- 
rameter components at the wavevectors in the star of 
Q we use the elements {£, C4, Cf, Cf} (these give the 
star of Q, {Q 1 , Q 2 , Q3, Q4} respectively, as shown in 
Fig. 1). This then defines a superconducting order pa- 
rameter with four components which we define as A-p i = 
(A ri , Ql , A ri ,Q 2 , A ri ,Q 3 , A ri ,q 4 ). With these definitions, 
the symmetry properties of the order parameter are given 
as follows (Dr^g) defined in Table I): translation T, 

Ar„Q 3 - < Q T A, .q ( Af ijP . - e~^- T A? 
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TABLE I: Representative spin-singlet, spin-triplet, and spin 
density basis functions for the different irreducible represen- 
tations that have momentum Q 1 = (q, q, 0.5). 



the transformations Gq lead to 

C 4 : D Ti (C i )(A T ^Q 21 A rii g 3 , A r ^Q 4 , A Ti Ql ) 

v z ■ D ri (a z )(A ri . Ql , A ri! q 2 , Ar 4 ,Q 3 , A r ,,Q 4 ) 

C2r, : -Dr i (C 2 , ) )(A ri: Q 1 , A ri ,Q 4 , Ar i: Q 3 , A ri ,Q 2 ) 

cr c : Dr 4 (a- c )(A rj! Q 1 , A ri ,Q 4 , Ar 4 ,Q 3 , A ri ,Q 2 ) 



(1) 



Table Q] reveals that both singlet and triplet order pa- 
rameters belong to the same representation which im- 
plies that singlet and triplet superconductivity are mixed. 
This mixing is due to spin-orbit coupling which cannot 



be justifiably ignored in CeCoIns. Formally, this is a con- 
sequence of the fact that Gq does not contain inversion 
symmetry [bH ]. Previous studies have examined the de- 
velopment of sing let-triplet mixing in related situations 



17L [18| , often through Lifshitz invariants that 



appear in the Ginzburg Landau free energy when parity 
symmetry is broken. We have confirmed that our formal- 
ism yields the same results as through the use of Lifshits 
invariants. 



Free Energy and PDW solutions: We use a Ginzburg 
Landau theory to describe the PDW and d-wave order 
parameters. While this will not be reliable on a quan- 
titative level, it will allow us to correctly identify the 
properties of the PDW order and make robust experi- 



mental predictions. The PDW Ginzburg Landau free en- 
ergy density is constructed by imposing invariancc under 
the above symmetries (note that this free energy density 
is the same for all Ff), 



/ = «E ( \±r uQi ? + &(£« |A ri , Q ,| 2 ) 2 +/? 2 £ 1<J |A r „ Qi | 2 |Ar„ Q J 2 + /? 3 (|A ri , Ql | 2 |A r! , Q3 | 2 + |A ri , Q2 | 2 |Ar,,Q 4 | 2 ) 
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where D = — zV — 2eA, B — V x A, D\ corresponds to 
the (1,1,0), and D 2 to the (1,-1,0) direction. The free 
energy density for the d-wave order parameter is taken 
to be: 

f d = a d \A d \ 2 + /3 d \A d \ 4 + K \DA d \ 2 + K c \D z A d \ 2 (3) 

The coupling between these order parameters is given by 
(note that this coupling is the same for all Ti): 

fc = &iEJA rf | 2 |Ar ()Qi | 2 

+/3 c2 [A 2 (Ar 1 , Ql A r! ,Q 3 + A ri ,Q 2 A ri ,Q 4 )* + c.c] 

(4) 

This free energy is similar to one studied earlier in 
the context of PDW order in Ref. Jl2]. The "homo- 
geneous" phase in the absence of a magnetic field has 



five PDW states distinct by symmetry, if we ignore the 
d-wave phase (the phase factors (j>t,<f>2, and 4> 3 are not 
determined by the free energy): 



= (e*, 0,0,0) 

= (e*,e*,0,0) 

= (e*,0,e*,0) 

= (g*"^i g*"fe gi</>3 g*(0i+<Aa— <Aa)^ 

= (e*^ 1 ie ic ^ 2 e^ 3 ig i W 1 +^' 3_< ^ 2 )) 



(5) 



This set is reduced to A[?^ and A^ in the presence 
of a d-wave order parameter (l2T |. Finally, a magnetic 
field in the basal plane along the (1, 1, 0)-direction favors 
the state Ap ; as it removes the degeneracy between the 
Q 1 and Q 2 wavevectors (the pairs Qi,Q 3 and Q 2 ,Q± 



3 



remain degenerate). As an example we assume that the 
field lies along (1, — 1, 0)-direction and yields the state 
(0,e*,0,e^ 4 ). Choosing cj> 2 = 4 = 0, the spatial de- 
pendence of the PDW order is given by A ri! g 2 cos(Q 2 • 
R) . In view of the coupling to the d-wave order param- 
eter the relative phase between A^ and Ar,,Q 2 can be 
either (jr) or ±7r/2 [ijj which are both permitted by 
the free energy. Generally, the combined PDW and d- 
wave superconductivity must then take one of two forms 
when vortices are ignored: + i| Ar ; ,2 cos(Q 2 -R) (time 
reversal violating phase) or A<2+|Ar,, 2 | cos(Q 2 -R) (time- 
reversal invariant phase) . 

Coupling to spin-density wave: The SDW order can 
be induced through the combined PDW and <i-wave su- 
perconductivity [3| . We assume here that the SDW is 
sufficiently weak so as to not alter the free energy sig- 
nificantly. The free energy density for the SDW order is 
then simply f S DW = a s S^ 2 St Q2 + f coupling with a s > 0. 
To determine f coupling, it is important to note that the 
SDW order breaks time reversal and can be generated 
by the PDW and d-wave order in two ways. The first 
is by coupling directly to the time-reversal symmetry 
violating phase and the second is by coupling to the 
applied magnetic field and the time-reversal invariant 
phase. This leads to two possible coupling terms that 
appear in f coupling, the first exists without the magnetic 
field, 

7iiS£ 3 {A2Ar 11<34 - A d A? liQ J + c.c. (6) 

and the second exists only in a finite magnetic field and 
is 

72 # 1 S§ 2 {ASA r4)Q4 + A d A* VijQ2 } + c.c. (7) 

where we have included Hi, the magnetic field along 
the (1, —1,0) direction. The experimental observation of 
a non-zero Sq 2 therefore leads to two possible types of 
PDW order. In the time-reversal broken phase, the PDW 
order must belong to the Ti representation. In the time 
reversal-invariant phase, the PDW order must belong to 
the L4 representation. This second possibility is most 
closely related to the 7r-triplet staggered phase that has 
been found within a simple microscopic description of 
CeCoIn 5 Q. Note that, in principle, both the represen- 
tations Ti and L4 will appear simultaneously. However, 
it is reasonable to expect that one of the two repre- 
sentations will give rise to the dominant order parameter. 

Role of Vortices: Now we examine the influence of the 
vortices on the induced SDW order, in the mixed phase 
in a magnetic field. Prior to turning to the detailed anal- 
ysis, we present the two main results here: (i) The vor- 
tex cores of the two PDW degrees of freedom A-p i ,Q 2 and 
Ar i! q 4 can lie at different positions and also need not co- 
incide with the d-w&ve vortex cores. We find that there 
exist stable phases where this happens. These phases are 



defined by the relative displacements r, of the PDW vor- 
tex cores from the d-wave vortex cores. In such phases, 
the d-wave vortex cores exhibit CDW order, (ii) The 
SDW order leads to Bragg peaks that are determined by 
the reciprocal lattice vectors of the vortex lattice and the 
displacements Ti (see Eq. 14). 

For a detailed derivation of the above results, we carry 
out the simplest realistic analysis. We assume that 
the correlation length of the spin-density order is much 
smaller than the coherence length of the superconducting 
order. We take Eq.[7]as the term driving the SDW order 
(the same arguments can be applied if Eq. [6] is used). 
From this we obtain 

S Z Q2 (R) = ^[A d (RrA ri , Q2 (R)+A d (R)A* i!Qi (R)}. 

(8) 

The spatial dependence of the PDW and d-wave or- 
der parameter can now be determined in the high-field 
limit for which the field H may be considered uniform. 
From Eq. [31 one finds that the d-wave component yields 
an Abrikosov vortex lattice. Using z to represent the 
(0,0,1)- and x the (1, 1, 0)-direction, the vortex lattice 
solution can be given by 

A d (x, ~z) = A d0 £ ^Uln-i/We-V-**)'/* (9) 

n 

where x — x/e, z = ez, the vortex lattice in the co- 
ordinates x, z has the basis vectors a — (a, 0) and 
b = 6(cosa,sina) 0, c„ = e ™pn 2 e -i*p(n+i) } q = 
2ir/a, z n = bsma(n + 1/2), p = (6/a)cosa, and 
e = [(k — K c )/n] 1 ' i . The parameter e scales lengths in 
the x and z directions to take the anisotropy into ac- 
count. This solution is an n — eigenstate of the op- 
erator D = D 2 X + D 2 Z — (— iV - 2eA) 2 with eigenval- 
ues (2n+ l)/l' 2 and I 2 = $o/(2ttH) ( n = 0,1,2,.. is 

the Landau level (LL) index). The macroscopic degen- 

- 2 

eracy of the eigenstates of D is exploited to create the 

Abrikosov vortex lattice solutions and, at the same time, 

plays a central part in constructing degenerate solutions 

for the displaced vortex lattice (4> n ) characterized by a 

vector r: (j) n (r +t) = e~ lTyX '<j> n (r + r) with <j) n {r) being 

a vortex lattice solution in LL n. The states <p n and <j) n 

~ 1 

are degenerate eigenstates of the operator D . 

In order to determine the PDW vortex structure it suf- 
fices to consider the linear equation for the PDW order 
parameter, which is found by keeping both Eqs. [3] and 
HI and by setting /?, = in Eq. [3J As a technical sim- 
plification, we set (ki — IK2D/K3 = k/(k + n c ) to ensure 

that the <i-wave order and the PDW order share the same 
~ 2 

D operator and hence have the same eigenstates (results 
without this simplification will be given elsewhere) . Min- 
imization of the free energy yields the following for the 
two degrees of freedom in the PDW order: 

nA r4 ,Q 2 = -/? cl |A rf | 2 A r4iQ2 - /3 c2 A* A*^ 



4 



nA 



r 4 ,Q4 



-/3 cl |A d | 2 A r4 , Q4 -/3 c2 A2A* 



,0a 



(10) 



with II = (a + \J («i — k 2 )(k! + k 3 )£> ) . To solve these 
equations, we expand the PDW order in eigenstates of the 

~ 2 

D operator. At sufficiently high fields, the PDW order 
will lie predominantly in the n = eigenstate for both 
Ar 4 ,Q 2 and Ar 4 ,Q 4 , and we ignore the smaller higher 
n contributions here. As mentioned above, these solu- 
tions are degenerate, implying the use of two displace- 
ment vectors t 2 and T4. At the second order transition 
where the PDW order appears, the vortex lattice struc- 
ture is determined entirely by the <i-wave order param- 
eter, so the only undetermined parameters are r 2 and 
T4. Solving the resulting linear equation yields the re- 
sult that the optimal PDW state is found by minimizing 
PciPa{t2) - |/3 C 2/3(t 2 , 1-4)1 with respect to t 2 and r 4 
where 



Mr) = ' ? e 

G 

/3(t 2 ,t 4 )= 

G 



iGT 



(11) 

(12) 



where G are the reciprocal lattice vectors of the vortex 



lattice, G = G- 



2-kB 



xt 2 and f3 = unless t 2 + t 4 = T, 



where T is a vortex lattice translation vector. For /3 cl < 

it follows immediately that r 2 — T4 — while the 
solution for /3 c i > requires a numerical minimization to 
determine t 2 . The resulting configurations are shown in 
Fig. 2, assuming that the <i-wave order forms a hexagonal 
vortex lattice. 

Here we provide a description of the phases in Fig. 2. 
The phase diagram depends upon r = |/3 C 2 //3ci | and in all 
the phases we can choose T4 = — t 2 . We find that there 
are four phases: in Phase 1 (0 < r < 0.07), 1-4 = 7(0 + 6) 
and 7 evolves continuously from 1/3 to 1/2; in Phase 2 
(0.07 < r < 0.31) 7 stays fixed at 1/2 (Fig. 2 shows 
r 4 = a/2 which is equivalent solution to x 4 = (a + b)/2); 
in Phase 3 (0.31 < r < 0.5), T4 = 7 2 a where 7 2 changes 
continuously from 1/2 to 0; finally in Phase 4 (r > 0.5), 
T4 = 0. The arguments of Ref. [l2j imply that in Phases 

1 through 3, the <i-wave vortex cores have charge density 
wave order at twice the PDW wave-vectors. 

The solution of the vortex lattice problem for the PDW 
order allows the SDW order to be determined which is 
particularly important as neutron scattering measures 
the Fourier transform of S Z (R). Eq.(fT2|) yields the in- 
triguing result that the SDW order will exhibit Bragg 
peaks at k positions that depend upon r 2 and T4: 



k = Q 2 + G 



2nB 



(13) 



where G is a reciprocal lattice vector of the vortex lattice 
and r is either t 2 or T4. Consequently, the relative po- 
sition of the vortex cores of the PDW and d-wave order 
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FIG. 2: Possible vortex configurations for the PDW order. 
The yellow circles give the zeroes of the d-wave order param- 
eter, the blue diamonds give the positions of the zeroes of 
Ar 4 ,Q 4 , and the red triangles give the positions of the zeroes 
of A r 4 ,Q 2 . 



can be retrieved from the position of the Bragg peaks in 
the SDW order. 

Conclusions: We have developed a phenomenological 
theory for the Q-phase of CeCoIn 5 to identify the possible 
symmetries for the PDW order. This theory is used to 
determine phases in which the PDW and d-wave vortex 
lattice are relatively displaced, leading to CDW order in 
the d-wave vortex cores. Interestingly, these structures 
can be probed by the position of the SDW Bragg peaks. 
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